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J>^ Abstract 

Motivated by their effect on the gravitational wave signal emitted by cosmic strings, 
we study the dynamics of kinks on strings of different tensions meeting at junctions. The 
propagation of a kink through a Y-junction leads to the formation of three 'daughter' 
kinks. Assuming a uniform distribution of the incoming wave vectors at the junction, 
we find there is a significant region of configuration space in which the sharpness of at 
least one of the daughter kinks is enhanced relative to the sharpness of the initial kink. 
Q_i' For closed loops with junctions we show this leads to an exponential growth in time of 

very sharp kinks. Using numerical simulations of realistic, evolving cosmic string loops 
with junctions to calculate the distribution of kink amplitudes as a function of time, we 
show that loops of this kind typically develop several orders of magnitude of very sharp 
kinks before the two junctions collide. This collision, or other effects such as gravitational 
VO \ backreaction, may end the proliferation. 

^ : 1 Introduction 

o. 

In recent years there has been a revival in the study of cosmic strings and cosmic superstrings, 
motivated largely by the realization that cosmic superstrings can form as a by-product of brane 
inflation (see e.g. [U El El IH E] for reviews). Of crucial importance is the fact that their 
observational signatures may provide a unique window on string theory. Hence much work has 
been dedicated to studying the evolution of cosmic superstring networks (see e.g. [SI El El CDl 
[IT] ) and to determining e.g. their gravitational wave signatures [121 [TBI Q31 [T3], which may be 
detectable by future experiments such as LISA. 

The properties of cosmic superstrings differ in at least three ways from those of standard 
cosmic strings, whose evolution and observational signatures have been studied in depth for 
over 30 years. First, a network of cosmic superstrings contains different types of strings, each 
with a different tension: fundamental F-strings; Dl-branes or D-strings; and (p, g)-strings which 
are a bound state of p F-strings and q D-strings. Second, a network of cosmic superstrings is 
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thought to contain numerous Y-junctions, namely points at which F- and D-strings meet to form 
the bound state (p, q) -string. Finally, whereas the inter-commutation probability for standard 
cosmic strings is P = 1, this is much reduced for cosmic superstrings [TB]. Indeed, for the 
collision of two F-strings, 10 _3 <Pp<l, whereas for the collision of two D-strings, 0.1<Pd<1. 
When strings of different types (such as an F-string and a D-string) collide, they cannot inter- 
commute due to flux conservation. In certain cases [UJ HE] they form a bound state string with 
two corresponding Y-junctions. 

Earlier studies on the detectability of standard cosmic strings through their gravitational 
wave emission have recently been partly generalized to cosmic superstrings [121 d31 EH]- In 
these papers cosmic superstrings were modeled as usual cosmic strings, but with a reduced 
inter-commutation probability P ^ 1. However, until now, the effects of junctions and bound 
states, as well as the implications of having strings with different tensions, have not been taken 
into account. 

In this and in a companion paper [TJJ], we study the effect of junctions on the gravitational 
wave burst emission from cosmic string networkgj. It is well-known that gravitational radiation 
from cosmic string loops is dominated by the lowest frequency modes (which are a multiple of 
the fundamental frequency for loops without junctions) [20] • Superimposed on the stochastic 
background of gravitational waves they generate are high frequency bursts emitted at cusps and 
kinks [21]. Kinks in particular radiate as they propagate along a string and through a junction, 
and when they interact with other kinks [22] • Cusps on the other hand are punctual in time, 
but generate bursts with a higher amplitude. 

For cosmic strings without junctions, cusps provide the dominant contribution to the GW 
burst signal. In [23] cusps have been argued to be a generic feature on strings with junctions 
as well. More recently, however, it has been shown [21] that since cosmic superstrings evolve in 
a higher dimensional space-time, cusps may be very rare events and furthermore those cusps 
which are formed are rounded off, hereby significantly reducing their GW burst signal. 

For this and other reasons we focus on kinks in this paper. On a loop with no junctions, 
the number of kinks is constant, fixed by the initial configuration of that loop. However, for a 
loop containing junctions this is no longer the case. Such loops evolve non-periodically in time 
and, as we will see, the number of kinks on them increases rapidly: kinks proliferate. Here we 
calculate how the number of large amplitude (that is, very sharp) kinks proliferates as the loop 
evolves, since these will turn out to dominate the GW signal from kinks on loops with junctions 

The paper is organized as follows. In Sections [2] and [31 we focus on the interaction between 
a kink and a single junction. (That is, we do not as yet consider closed loops - - this is 
done in Sections H] and [SJ ) In Section [2] we use the dynamical equations of motion for strings 
with junctions derived in [17J to show that when a kink propagates through a Y-j unction, it 
leads to the formation of three 'daughter' kinks (one reflected, and two transmitted). For a 
specific junction configuration for which the whole evolution can be solved analytically, we 
show explicitly that the amplitude of the daughter kinks may be larger than that of the original 
'parent' kinkjj In section [3] we generalize this discussion by considering arbitrary junction 
configurations. More specifically, we take a uniform distribution of incoming waves at the 
junction and of incoming kink amplitudes and show that, in a significant region of configuration 
space, the amplitude of at least one of the daughter kinks is larger than that of the parent kink. 



^^Note that spontaneous symmetry breaking phase transitions can also lead to the formation of cosmic strings 
with junctions, and hence our analysis also applies to such strings. 

2 Kink amplitude, defined in section [H is synonymous to the kink sharpness - a nomenclature often used 
elsewhere. 



In Sections H] and El we study the evolution of the number of sharp kinks on closed loops with 
junctions. In Section HI we consider a simplified model of a loop with junctions that does not 
take into account the complicated dynamics of the loop itself, and we argue that the number 
of large amplitude kinks increases exponentially with time. Finally, the dynamics of the loop is 
taken into account in Section El where we show that this dynamics generally further enhances 
the proliferation of sharp kinks. Our conclusions are presented in Section El 

2 Propagation of a kink through a junction 

In this section, following |17j . we first review the equations of motion for three semi-infinite 
Nambu-Goto strings of tensions fii, \ii and /X3 which meet at a Y-junction. A consequence of 
these equations is that the propagation of a kink through a Y-junction results in the production 
of three 'daughter' kinks; a reflected kink as well as two transmitted kinks. We also define kink 
amplitude (or sharpness), and analytically study - in the simplest case of an initially static 
junction - the amplitude of the daughter kinks, showing that in some cases these can be 
amplified relative to the incoming kink. 

2.1 Description of the system 

We work in flat spacetime with signature ( — h ++), and use the standard conformal-temporal 
gauge so that each string is described by its spatial coordinates Xj(cr, t), where t coincides with 
Lorentz time and the subscript j labels the different strings. The gauge constraints can then 
be written as (with ' and ' standing for derivatives with respect to a and t respectively) 

x ;.-x, = (i) 

x? + x? = 1. (2) 

The action describing the combined system of three semi-infinite strings meeting at a junction 
has been analysed in [T7j. Away from the junction the wave-like equation of motion for each 
string has solution 

x jM) = 2( a i( M )+ b i(* ; )) ( 3 ) 

where 

u = a + 1, v = a — t, (4) 

and where a'? = b' =1 in order to satisfy the gauge constraints. Each string is bounded by 
the junction located at X(t) = Xj(sj(t),t). One can therefore let a take values in the interval 

] - OO, Sj(t)}. 

As explained in [T7], the initial conditions for Xj and x'- at t — determine slj(u) and bj(v) 
for u ^ Sj(0) and v ^ Sj(0). A set of coupled differential equations describing the physics at 
the junction then enables one to determine the evolution of Sj(t) as well as the outgoing waves 
a' (sj(t) + t) in terms of the incoming waves b'Asj(t) — t) (which are determined by the initial 
conditions since Sj ^ 1). In particular, let 

v\ = /J-2 + A*3 - Mi > (5) 

Mi = /il- (fjQ - /i 3 ) 2 = ^3 ^ 
c x (t) = b' 2 (s 2 (t)-t)-b' 3 (s 3 (t)-t), (6) 



and circular permutations, as well as 

jJ, = (ii + /i 2 + /i3- 

Then the equations of motion imply energy conservation at the junction 



(7) 



One also has 



Hj{l 



/Mi + H2S2 + H3h = 0. 



MJl-Cj 



/' 



and 



1 ;.( Sj (t) + 



Mi(l - ci) + M 2 (l - c 2 ) + M 3 (l - c 3 ) 



f 1 fc=i 



(9) 



(10) 



which determines the outgoing wave on string j in terms of the inward moving waves. The last 
term in (fTUj) is proportional to the velocity of the junction: 



X 



■-EMi-^Kk(t)-t). 



f \\\ 



From (HOP it follows that if one of the strings has a kink propagating towards the junction, 
i.e. one of the functions b'(t>) has a discontinuity, then all a'-(tt) acquire a discontinuity when 
the kink reaches the junction. The presence of the junction therefore increases the number of 
kinks in the system from 1 to 3. We refer to the three newly formed kinks as 'daughter kinks'. 

Furthermore, eq. (1101) also implies that for essentially all initial conditions fixed at t — 
say, the subsequent evolution always generates kinks. Indeed, from Eq. (fTUj) and for t = + , 
&'j(sj(0) + ) on string j is determined by the weighted sum of all b' k (sk(0)~). The latter are fixed 
by the arbitrary initial conditions and independent of the &' k (sk(0)~). Their sum therefore 
yields a vector a'(sj(0) + ) which is generally different from a'-(sj(0)~). Therefore a' will be 
discontinuous at u = Sj(0), leading to a kink t = + . Hence the presence of junctions essentially 
implies the existence of kinks, which therefore will not need to be introduced by hand in our 
simulations below. 



2.2 Amplitude and transmission coefficients 

Kinks on cosmic strings are sources of gravitational wave bursts: a kink emits bursts throughout 
its propagation on the string [2T], and also when it encounters another kink or when it crosses 
a junction [22]. However, the GW signal emanating from kinks on a network of strings is 
determined not only by their number, but also by their amplitudes. The kink amplitude can 
be defined as follows. 

Consider a kink moving towards a junction (i.e. a discontinuity of h'Av) at u # ). The ampli- 
tude of all GW bursts associated with kinks is proportional to the components of [2TJ 



b'l 



mv; 



l-n.b£(v+) l-n-tyfo-) 



(12) 



3 The exception is when, given the h'j(sj(0) ), the a^(sj(0) ) are chosen to be given by the RHS of eq. |T 
A particular case of this are totally static initial conditions (xj = = X = Sj) considered in [25j . 



Both denominators in ( 1121) are generically comparable and of order one, because the direction 
of emission n is a priori uncorrelated with b'(u+) and b'(w~), and therefore (1121) is of order 
l|ko-(^*~) — k'"( v «T)ll- We define the kink amplitude or sharpness by 

A[h' 3 ] = i||ty(t;+) - b;.(0|| = | sin(0/2)| (13) 

where 9 is the angle between h'Avf). The factor of 1/2 is a normalisation factor so that 
< A[h'A < 1. We adopt an analogous definition of the sharpness of an outward moving kink 
on string j characterized by a discontinuity in a'-(u). 

At first sight one might expect the daughter kinks to have smaller amplitudes than the 
incoming kink. However, this is not the case. In particular energy conservation at the junction, 
eq. (El), does not constrain the amplitude of the transmitted kinksj. Indeed the derivatives in 
(JSJ) suffer a discontinuity when a kink hits the junction, and these therefore undergo a sudden 
jump. 

It will also be useful to define the following transmission coefficients Cj. For instance, 
consider a kink that propagates towards the junction on string 1 (i.e. h' ± (v) has a discontinuity 
at vu) and reaches the junction at £*, with vu = s\(U) —t*. Then each function a^-(-u) acquires 
a discontinuity at u^ = Sj(t*) + £*. We define 

C - 431 (14) 

'~ A\h\]' (14) 

where the amplitudes are given by eq. ffT31) . 

2.3 Example: static junction 

To conclude this section, we illustrate the production of daughter kinks with an example that 
can be worked out analytically. 

The strings of tension fi\ and fi2 = /^3 are taken to be in the (x, y) plane at all times, and 
the initial configuration considered is shown in the left-hand panel of Figure 1. String 1 lies 
along the x-axis, and strings 2 and 3 subtend an angle ip w.r.t the |/-axis, chosen such that the 
junction is initially static. There are two kinks on string 1, both of which propagate towards 
the junction (they are discontinuities in h[). Once both kinks have propagated through the 
junction, the angles between the strings at the junction is again the same as initially and hence 
the junction is again static. However, all three strings now have outward moving kinks on them 
as shown in the right hand panel of Figure 1. The amplitude of the outgoing kinks can be 
calculated as follows. 

Let 

ki = (-1,0) 

k? = -(cos#,sin#) 

k 2 = (e,-6) 

k 3 = (e,8) (15) 

where 9 G [— n, n] is a free parameter and we have defined 

0<e=^-<l, 5 = VT^^. (16) 

2/i2 



4 Thc energy along the string is proportional to a. Hence no energy is associated to the kink itself which is 
described by a pointlikc discontinuity. 



Hence the angle ip = arctan(e/5). Then the initial condition corresponding to Fig. 1 is given 
by 











k! 


if - L ^ f ^ 




ai(«) = 


= k x 


if u < 


KM = 


k? 


if - L - £ < v < 


-L 


a 2 (w) = 
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= k 3 


if u ^ 


b' 2 („) = 
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From 



it now follows that for string 1 and for v G [-L — £, —L], 

;(0/2)(cos(0/2),sin(0/2)" 



cos 



*i = sin(0/2)(sin(0/2),-cos(0/2)) . 



(17) 



Thus the physical angle to the x-axis made by the segment of string 1 between the two kinks 
is # g eo = 0/2. It moves with velocity sin(0/2) towards the junction. 





'M3 = M2 



Figure 1: Propagation of two kinks through an initially static junction. Initial (left) and final 
(right) configuration of the system. 

From the equations of motion (ioT)-(TTUl) one can analytically determine the evolution of the 
system, since the b' are piecewise constant functions. The evolution consists of three distinct 
phases: 

• te [o,L]. 

The solutions of Cj in (JBJ), substituted into (jHD, yield Sj = on each string. Thus ffTTT) 
implies the junction is indeed static: X = [IT]. The outgoing waves are determined by 
(1101) : a'Au) = kj for u G [0,L]. Hence during this phase, the segment between the two kinks 
propagates towards the junction with velocity | sin(0/2)| whilst the rest of the system remains 
static. This phase ends when t = L, because then si(t = L) — L = v = —L so the scalar 
products Cj(t) (which depend on b' 1 (si(t) — £)) change. 

• t G [L,Q. 

During this phase the segment between the two kinks crosses the junction and b^ = Iq. 
The time i* is the solution to Si(t*) — t* = —L — £. Beyond this time, b[ changes back to ki 



again. From ([HI) we find that the Sj are constant and equal to 



e (cos9 — 1) 

Sl = ecos# + (l + e + e 2 ) ' (18) 

(l-cos9)e 2 -5sm9(l + e) . . 

s 2 = /i , /i , T^n ' ( 19 ) 



ecosfl + (l + e + e 2 ) ' 
(l-cos#)e 2 -5sin#(l + 





(1 


e cos6> + (1 +e + e 2 ) 
-cos6>)e 2 + 5 sin 9(1 + 


e) 



e cosfl + (1 + e + e^J 

Thus Sj(t) = Sj(t — L), so that t* = L + jzj~- The functions a'- can then be obtained from (1101) : 
for instance 



a ' l(B) = l+^ + 2 t coso ( " (1 + f2)COS " " 2f ^ 2sin ») " e 



L,L+ 1 + *' 



1-s. 

Finally eq. (1111) determines the constant velocity of the junction during this time interval; 



(21) 



X = — I ^(costf - 1), (1 + esinfl)) . (22) 

1 + e + e z + e cos v 7 



t > t m . 



The configuration at the junction is now exactly the same as during the first phase. That 
is, Sj = on each string and b!j(u) = kj for u > L + jz$-- Now, however, there is a segment 
between two kinks on each string and it propagates away from the junction as depicted in the 
right panel of Fig. [TJ 

Calculation of the transmission coefficients (j!4p when the first kink encounters the junction 
yields 

Ci(0,e) = ; (1 ~ e) = <1, (23) 

V ; V2ecos^ + e 2 + l V ; 

C 2 (M = 2£ , (24) 

^(1 + e + 2e 2 ) + e(l - e) cos^ - (1 + e) v / T^ 2 sin^ 

C 3 (^e) = C 2 (-0,e). (25) 

In Fig. El we plot the transmission coefficientqj for different values of the string tensions in the 
allowed range ^ e = /ii/2/i2 ^ 1. One sees that, in this particular example, the reflected kink 
on string 1 always has a smaller amplitude than the incoming kink, though the reduction in 
amplitude is generally rather weak for incoming kinks on the lightest string. 

By contrast, the transmitted kinks can be amplified. When all tensions are equal (e = 1/2), 
this occurs for a rather broad set of (static) junction configurations. However for e — y (that 
is, when strings 2 and 3 are heavy compared to string 1), we find C% > 1 only in a limited range 
of 9. In this regime, it is in fact straightforward to understand the position of the peak. From 
( 1101) it follows that the amplitude A\a!A - and hence the transmission coefficient Cj - is large 
when after the kink has crossed the junction, the corresponding Sj — > — 1. Eq. (ITU1) predicts 
that, for small e, this occurs on string 2 when sin^ = 1 or 9 = 7r/2, as is indeed the case in 
Fig. El Away from its sharp peak, Ci ~ e for small e, as can be seen from eq. ( 1241) . Finally, for 
e ~ 1, Ci is always close to 1 even though slight amplification can occur in a broad range of 9. 



3 We only plot C\ and C2, since C^(6) = C2(—0) because of the symmetry of the initial configuration. 
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Figure 2: Transmission coefficients C\ and C 2 as functions of 9 for three different ratios of 
tensions \i\ and \i 2 - The tension /i 3 of the third string is taken to be equal to [i 2 in all three 
cases. 

We should stress that the above discussion is limited to the specific junction configuration 
of fig. [TJ We now turn to arbitrary junction configurations. 

3 Distributions of transmission coefficients 

In this section we aim to gain intuition on the distributions of transmission coefficients char- 
acterizing the propagation of kinks through junctions in a more generic context. Namely we 
study the statistical properties of the Cj using an underlying uniform distribution of junction 
configurations. We would like to answer the following question: is amplification frequent or 
not? 

To specify the configuration of the junction just before the arrival of the incoming kink as 
well as the amplitude of the kink one needs four incoming unit vectors. We consider a kink 
moving towards the junction on string 1, specified by a discontinuity in b[(v) at vu- Let t* 
be the time when the kink reaches the junction. Then the amplitudes of the transmitted kinks 
depend on b'^uj^), b'^uj^), b' 2 (s2(t*) — £*) and bg^^*) — £*), for which we will use the more 
concise notation b' x + , b'{~, b' 2 and b' 3 . Eq. ffTU]) then yields 



„/± 



a? 



± u'± 



ffbf - Qt 2 K - Qt 3 b' 3 , 
P 2 b 2 — Q 2i3 b 3 - 
V3,ibi 



PtK 



Qtibi 
- Q 32 b r 



(26) 

(27) 
(28) 



where 



>± 



and the sf are given by eq. 




Q 



± 



Using eqs. 




(29) 



the transmission coefficients Cj can be 



obtained from 014|) . Our aim here is to determine their probability distributions. Here we focus 
on the case of equal tensions. Unequal tensions are studied in Appendix A. 

The configuration of the junction just before the arrival of the kink is specified by the 
unit vectors b'f , b' 2 and b' 3 . We assume for now that these are independent, with uniform 
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distributions on the unit sphere. We also assume a flat distribution in the incoming kink, 
namely a uniform distribution on the unit sphere for b'^ (and in particular independence of 
b'f 1 "" of the other unit vectors jj. 

Given these assumptions, we numerically calculate the distributions of the transmission 
coefficient by drawing a large number (typically JV = 3x 10 7 ) of random configurations at the 
junction in order to estimate the various statistical quantities of interest. 



3.1 Marginal Distributions p(Cj) 

We are in the first place interested in the joint distribution p(Ci, C2, C3) from which one could 
determine, for example, the probability for several kinks to be amplified at the same time. 
However, as a warmup, we show in Figure [3] the marginal distributions p{Cj), for the case of 
equal tensions fij = 1. These distributions p{Cj) have a large tail where Cj > 1, indicating 
that both the reflected and the transmitted kinks can be amplified in a significant part of 
configuration space. Indeed, P{C X > 1) = 0.11, P(C 2 > 1) = 0.19 = P(C 3 > 1). 




Figure 3: Marginal distributions p{Cj) of the transmission coefficients Cj for equal tensions 
fij = 1 and a uniform distribution of junction configurations. The average values are (Ci) = 
0.49, (C 2 ) = 0.72, (C 3 > = 0.72. 



However, as expected from the static example of Section 12.3} the marginal distributions 
depend significantly on the ratios of string tensions. In Appendix A we discuss the distributions 
for various other sets of tensions, including several limiting cases where, using the smallness of 
some of the coefficients (12"9"1) . analytic arguments can explain certain features of the distributions. 



3.2 Simultaneous amplification of transmission coefficients 

In order to understand how the total number of large amplitude kinks on strings containing 
Y-junctions changes in time, one really needs the joint distribution p(Ci, C 2 , C 3 ) which contains 



6 The consequences of the first hypothesis have been studied in ]2<g\ where, for example, the probability 
distribution of the s~ was calculated. We will see in the next Sections that the second hypothesis does not hold 
for loops. It will turn out that most kinks on loops have a small amplitude. The flat distributions we use here 
should be interpreted as a working hypothesis adopted for now in the absence of a concrete dynamical model. 
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information about the correlations between the different transmission coefficients. 

The simplest question to ask involving correlations between the Cj is the following: when 
a kink reaches a junction, what is the probability that at least one of the three daughter kinks 
is amplified (so at least one of the Cj > 1)? This probability may well be significantly larger 
than that suggested by the tails of the individual daughter kink distributions discussed above. 
In the case of equal tensions fa — 1 corresponding to Fig. [3] we find 



P(at least one amplification) = 0.43 for fij = 1. 



(30) 



Amplification is therefore not such a rare event, and hence we have a first hint that the number 
of large amplitude kinks may grow significantly in a system in which the total number of kinks 
increases due to the presence of Y-junctions. 

In a similar way, the probability of having at least two simultaneous amplifications is 



P(at least two amplifications) = 0.07 for \ij = 1. 



(31) 



Again those probabilities are not negligible and such events can contribute importantly to the 
enhancement of the number of large amplitude kinks in an interconnected network^]. In the table 
below, we summarize the amplification probabilities for the different sets of tensions considered 
here and in Appendix A. 



tensions 


P{Ci > 1) 


P(C 2 > 1) 


P(C 3 > 1) 


P(at least 1 amp) 


P(at least 2 amp) 


fl! = 0.1, jU 2 = 1, M3 = 1 


0.12 


0.01 


0.01 


0.12 


0.004 


Ml = 1) M2 = 0.1, ^3 = 1 


0.01 


0.19 


0.19 


0.65 


0.006 


/Ul = 1.9, H2 = 1, A*3 = 1 


0.12 


0.47 


0.47 


0.93 


0.11 


Ml = 1, M2 = 1-9, M3 = 1 


0.01 


0.49 


0.01 


0.19 


0.004 


Ml = 1, M2 = 1, M3 = 1 


0.12 


0.19 


0.19 


0.43 


0.07 


/ii = 1, /i 2 = 1.2, ^ 3 = 1.4 


0.10 


0.11 


0.15 


0.31 


0.65 



3.3 Joint distribution p(C\, C2, C3) 

Even though there is a significant region of configuration space in which the amplitude of at 
least one of the daughter kinks is enhanced relative to the amplitude of the initial kink, it 
is important to know whether the amplitude of the remaining kinks is typically significantly 
reduced in such events. Fig. H] shows 2D slices through the joint distribution p(C\, C2, C3) for 
four different values of C±, again with fa = 1. (Note that the grey scale differs from one panel 
to the next.) The horizontal and vertical axes label C2 and C3 respectively. 
The joint probability distribution exhibits the following important features: 

• The top left-hand panel shows that there is a sharp peak where all transmission coefficients 
take small values, Cj <C 1. 

• The remaining three panels are for values of C\ > 0.2. They have a clear concentration of 
events on a arc-shaped line, as well as as events off that line at values of C2 ~ C3. These 
latter events saturate around C2 ~ C3 < 1 for C\ > 1. 

• The arc-shaped lines finish on two "amplification tails" corresponding to values of C2 or 
C3 larger than 1. There are no events for which Ci > 1 and C3 > 1 simultaneously. 



7 Numerically we have found no example of a simultaneous triple amplification (no matter how large the 
number of random configurations generated). We have checked that P(C\ > a, C2 > a, C3 > a) is non-zero 
any value of a strictly smaller than 1 (increasing the number of configurations always ends up yielding such an 
event) but this vanishies to zero when a = 1. 
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Figure 4: Slices of the joint distribution p(C\, C2, C3) for four different values of C\ and fij = 1. 
The horizontal and vertical axes label Ci and C3 respectively. The values of C\ on the different 
slices are C\ = (top left), C\ = 0.2 (top right), C\ — 1 (bottom left) and C\ = 1.5 (bottom 
right). Note that the grey scale differs from slice to slice. 

Based on the above properties of the joint distribution one can conclude that the bulk of 
the transmission events falls in one of the following categories; 

• All daughter kinks have drastically reduced amplitudes 

• All daughter kinks have slightly reduced amplitudes 

• The amplitude of the reflected kink is significantly reduced, one transmitted kink is am- 
plified and one is slightly reduced 

• The reflected kink is amplified and the amplitude of the transmitted kinks is comparable 
to the amplitude of the incoming kink 

We note in particular that the amplification of one of the daughter kinks does not imply 
that the amplitudes of the remaining daughter kinks are small. In fact, except for the first kind 
of transmission — which is not considerably more frequent than the others, see Fig. [3] — the 
amplitude of the scattered kinks is never strongly suppressed. This indicates that amplifications, 
combined with the rapid growth of the total number of kinks in a system with junctions, lead 
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to a large number of "large amplitude" kinks. The next section is devoted to a quantitative 
study of this phenomenon. 



4 Proliferation of large amplitude kinks on a loop with 
junctions 

We have seen that the sharpness of the daughter kinks is generally comparable and occasionally 
larger than that of the incoming kink into a junction. Motivated by GW physics, it is therefore 
interesting to study how the number of kinks with a large amplitude (of order 1) evolves on 
closed loops with junctions. Indeed, a sustained growth of the number of such 'large amplitude 
kinks' may well have an impact on the gravitational wave burst signal emanating from strings 
of this kind pS]. 

The total number of kinks on a closed loop with junctions increases exponentially in time: 
when an initial kink reaches a junction it gives rise to three daughter kinks, which in turn 
propagate towards another junction where they multiply again, and so forth. Even though 
the amplitude of most of the daughter and higher generation kinks is small, it is clear that 
the number of large amplitude kinks will also grow exponentially provided amplification occurs 
sufficiently frequently. Here we show this is indeed the case in a simple model of a loop with 
Y-junctions, based as before on an underlying uniform distribution of junction configurations 
at the time of arrival of the kink. The evolution of the loop, for a relatively general class of 
initial loop configurations, will be taken in account in section \5\ where we will see that the 
results obtained here remain largely valid. 



4.1 Setup 

The simplest example of a closed system with Y-junctions is a loop formed by three strings 
meeting at two junctions, a typical example of which is shown in FigO 



string 1 ,\i^ 
string 2,\i 2 
string 3,^i 3 




Figure 5: A typical loop formed by 3 strings and 2 junctions 
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Our loop model is based on the following two simplifying assumptions; 

1. all three strings joining the two junctions have essentially the same constant invariant 
length L; 

2. when a kink (whose amplitude is known) reaches a junction, the configuration of the latter 
is randomly drawn among those that yield the correct amplitude for the original kink. 

More explicitly we proceed as follows. First we specify the initial conditions at t — 0, namely 
K, the total number of left moving and right moving kinks on all three strings, as well as their 
amplitudes. (The initial positions of those kinks do not enter the subsequent analysis.) After 
an interval of time L, the first assumption ensures that all the initial kinks have reached a 
junction, but that none of their descendants have done so (recall that kinks propagate at the 
speed of light on the strings). Therefore, at t = L, the system contains 3K kinks, and at time 
t = nL with n an integer, the total number of kinks in the systenjj is 3 n K. We then evaluate 

Qj ( n ) = number of kinks on string j of amplitude > A 

as a function of the generation n, or equivalently time t = nL. 

We do so as follows. Consider a given (say inward moving) kink of amplitude Ai on string 
1. Following assumption 2 we determine the amplitude of the daughter kinks by randomly 
drawing junction configurations, similarly to Section^ but now we restrict the choice to those 
configurations for which ||t>^ — b'^H = Ai. (In practice, we first draw b'j - , b' 2 and b' 3 with 
uniform probabilities on the unit sphere and then draw b'/" with uniform probability on the 
intersection of the unit sphere and the sphere of center h'{~ and radius Ai, i.e. on a circle of 
radius Ai.) In other words, instead of using the probability distribution p(C\, C 2 , C 3 ) to draw 
the transmission coefficients, we use an improved version p^.(Ci, C 2 , C 3 ) that takes into account 
the amplitude of the incoming kink. 

To summarize, for each initial kink, our model amounts to building a tree of kinks. On this 
tree each node represents a kink, and contains the value of its amplitude as well as the number 
of the string on which it arrived. It has three daughter nodes whose values are drawn randomly 
according to the rules explained above. At the nth generation, the total number of kinks is 
3 n . Clearly, as there is no interaction between kinks, for K initial kinks there are K trees 
which evolve completely independently. Hence the statistical properties of systems originally 
containing several kinks can be trivially deduced from those by linearity. 

4.1.1 Numerics 

Computationally, as the number of generations increases it becomes expensive to store 3 n 
amplitudes. To keep the computation manageable we divide the amplitude interval [0, 1] into 
Nbi n bins, and only keep track of the number of kinks on each string with an amplitude in the 
different bins, at each generation. To draw the amplitudes of the subsequent generation we use 
the center value of the bin as the amplitude of the incoming kink. That is, for kinks in bin 
k > 1, we use as an amplitude (k — \)jr—- ( in our simulations, we used N bins = 100). Finally, 



8 In fact this is not exactly true if the constant invariant lengths of the strings are exactly equal. In this case, 
e.g. with K = 1 the second generation kinks reach the second junction simultaneously and therefore recombinc 
into a single kink per string instead of three. From then onwards, the total number of kinks in the system would 
remain constant and equal to three. We disregard this possibility and assume that, because of variations in the 
length of the strings when its evolutions is taken in account, the three kinks arrive at different times at the 
junctions giving rise to nine kinks. 
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we set the amplitude of kinks in the first bin to zero: we expect that after a few generations 
most of the kinks that lie in this bin actually have an amplitude smaller than 1/iVW by several 
orders of magnitude. Of course, this procedure can lead to an underestimation of the number of 
large amplitude kinks (because a few kinks from this bin must in reality be reamplified to yield 
large amplitude descendents) but it prevents many very small amplitude kinks from spuriously 
leading to large amplitude kinks. By doing this, we loose any information on the low amplitude 
part of the distribution (which we disregard in this paper). In our companion paper [IS], we 
will refine this numerical setup to demonstrate that the large amplitude part of the distribution 
dominates the gravitational wave signal. 



4.2 Results 

In FigE]we plot logQ/ (n) as a function of the generation n, for different sets of string tensions 
/ii,/i2 and /i3. (The particular choice of A = 1/4 will be justified in |19j.) The initial condition 
at the start of the simulation was a single right-moving kink K = 1 of maximal amplitude on 
each string. 
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Figure 6: The logarithm of the number of kinks on string 1 of amplitude A ^ 1/4, as a function 
of time measured in units of the constant string length L, for different sets of string tensions 
/ii,/i 2 and /i 3 . 

After some initial fluctuations at small n, which vary significantly from one realization to 
anothecl, one sees that the points corresponding to a given set of string tensions can be fitted 
to a straight line, whose slope is independent of the particular realization. Hence Q ± (n) grows 
exponentially with n. More generally we find 10 ! 



Qf(n) oc exp[7n] 



(32) 



9 This is because as long as the number of kinks is small, the random values of the Cj that are drawn crucially 
affect the distribution of kinks. 

10 Empirically we find the slope of logQ^(n) is independent of j for large n. 
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where the coefficient 7 depends on the tensions fij as well as on the amplitude A. One sees 
proliferation is most efficient when all tensions are equal. Gradually moving away from this 
case, the slope 7 = 7/ln(10) of the curves in Fig O decreases and approaches zero when one or 
more of the 17 given in (J5J) vanish. For a given ratio of tensions, we find 7 is approximately 
independent of A over most of the range of possible amplitudes. For equal tensions and A > 0.1 
shown in Fig E] one has 7 ~ 1/3. In the limit A — > 0, the slope sharply increases to 7 « In 3, 
since obviously Q°j{n) = 3 n . 

We therefore conclude that, at least in this simplified model of a loop with junctions, for 
a large range of string tensions, the amplification rate is sufficient to sustain an exponential 
growth of the number of large amplitude kinks. In appendix B we illustrate, with a toy ana- 
lytical model, the origin of this exponential growth. We must emphasize, however, that since 
amplification remains a rare event, the vast majority of kinks at sufficiently late times will have 
small amplitude. Indeed, the fraction of the total number of kinks that have a large amplitude 
tends to zero. The implications of these findings for the gravitational wave signal emitted by 
string loops of this kind will be studied elsewhere [T5] . 

5 Proliferation of large amplitude kinks on evolving loops 
with junctions 

The model of a loop with junctions discussed in sectionH]does not take in account the dynamics 
of the loop. In particular, it assumes the invariant length of the strings forming the loop is 
constant in time. Here we include the effects of the loop dynamics on the evolution of the 
number of large amplitude kinks, by numerically integrating the equations of motion ([§]) and 
( II Ul) . suitably modified to take account of the presence of the two junctions. 

This is more complicated than for periodic loops with no junctions, since one needs to keep 
track of the position of each junction, and extend the definition of the a' beyond their initial 
domain. Our simulations generalize those of [25] in which certain initially static and symmetric 
planar loop configurations were studied. Here we consider a rather general class of non-static 
initial conditions, and count the number as well as the amplitude of the kinks on the loop as a 
function of time. Note, however, that we do not take in account self-intersections between the 
strings in our simulations. This may be an important limitation of our model if the probability 
of intercommutation is large. (In this case, intercommutation between strings of the same type 
would lead to a loop being chopped off, whereas intercommutation between strings of different 
types might increase the number of junctions in the loop.) Our simulations end when the length 
of one of the strings connecting the two junctions shrinks to zero and the junctions collide. The 
outcome of such a collision is an open question, but it may well lead to the formation of two 
loops without junctions (assumed for e.g. in [27]): this would end the proliferation of kinks. 

5.1 Initial conditions 

Our initial condition consists of 3 segments of string of initial invariant length Lj(t = 0) 
(specified below) which join two junctions at positions 

X A = (0,0,-1) X fl = (0,0,1). (33) 

Motivated by the harmonic construction of standard cosmic string loops with no junctions 
I3"U| I3"TI I3"2l [33], the three initial string segments are taken to be three arcs of circles to 
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which we add a higher harmonic. Explicitly, 

/ - 

x,(a,t = 0)=n(a 3 ) 

V 



sin 



Hi 







COS 



11, 



l + bjSinlrijf- 



(34) 



/ 



where rrij,nj G N, and TZ(otj) is the rotation matrix about the z axis by an angle olj 

I cos(oj) — sin(aj) \ 
TI(ocj) = | sin(aj) cos(aj) ] . (35) 



srn(aj 




cos(aj 




If we set cij = bj = 0, then the segments are initially semi-circles. If Oj or bj are nonzero then 
the strings appear as perturbed semi-circles (at least for small values of a,j and bj), and the 
integers rrij and rij set the 'wiggliness' of the perturbation. The normalisation 



Hj = yj(l + a,) 2 + (a jmj ) 2 + (1 + bj) 2 + (bjUj)' 



(36) 



ensures that ||x'|| < 1 since x| + x'? = 1 from the gauge conditions (J2J). Note that with this 
definition of Hj in general sup ||x'-|| < 1 so that the initial condition is never static. The initial 
invariant length of the strings segments is Lj(0) = ttHj so that 



a E [-7cHj,0] at t = 0. 
In order to compare our results with Section 4, we will measure time in units of 



L = -(L 1 (t) + L 2 (t) + L 3 (t)), 



(37) 



(38) 



the total invariant length of the loop divided by the number of string segments. Note that this 
quantity remains constant throughout the evolution. 

The initial velocity of the strings Xj must satisfy the gauge condition Xj ■ x'- = 0. We choose 
it to be 



±j(a,t = 0) = Nj{a)K{ctj 



( sin (l?-) { l + b i sin [ n iir)) ~ h i n i cos (/^j) cos (™i]|) \ 



V - cos (/|-)( 1 + a i sm l m ^ 



anirii sin tt cos m,TT 



(39) 



u 3 " "3 
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with Nj(a) defined so that x 2 - = 1 — xf, and Vj is the component of velocity transverse to the 
plane of the string. Finally the functions &!j(z) and b'j(z) on the interval [—ttHj, 0] are obtained 
through 



a}(*) 
b'Az) 



x^a = z,t = 0) + Xj(cr = z,t = 0) 
x'(<7 = z,t — 0) — x.j(a = z, t = 0). 



(40) 
(41) 



As explained in Section El the cr parameter on each string takes values in [sA,j(t), SB,j(t)] at 
time t. Furthermore, integrating the equations of motion extends the definition of the functions 
b!Az) and b'Az) outside the interval [—nHj, 0] (using equations (TT0T) ) to z > for b!Az) and to 
z < —nHj for b'Az). Thus if the evolution is calculated up to a final time tf, then at the end 
of the simulation a'j(z) will be defined in the interval [—irHj, Zfj] where Zfj = tf + ss,j(tf)- 
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5.2 Proliferation of large amplitude kinks 



C-C -i - yJj i - t/ o - 



m,, n. 



and Wj, and 



This class of initial conditions has 6 parameters for each string: 
therefore enables one to probe a variety of initial configurations. We now evolve these strings 
and count the number of large amplitude kinks as a function of time. That is, we calculate 
A[a! 3 ]{z) = l\\a' j (z + ) - a^")|| (for left-moving kinks) and A[b$](*) = l\\h' j (z + ) - ty.(z")ll as 
functions of z. These functions are zero except at the position of a kink, where they reduce to 
the kink amplitude. 

As discussed in section [21 even though the initial configuration is infinitely smooth and 
appears to contain no kinks, this is not the case. While a' and b'- are continuous inside the 
interval ] —irHj, 0[ they have a discontinuity at z = ~^Hj for b'j(z) and z = for a!j(z). Indeed, 
as soon as the evolution starts, the equations of motion define the function a'-(z) for z > 0. 
In particular, b!Az = + ) depends only on the values of the h' e (z = 0") for £ = 1,2,3 and 
differ from a'Az — 0~). Our loop therefore initially contains 6 kinks: one left-moving and one 
right-moving on each string of the loop. 

The numerical integration of the equations of motion ends at time tf when the two junctions 
collide. The results are shown in figure [7] for the case of equal tension strings fij = 1. Here we 
have used the following set of parameters as initial conditions: string 1 (ol\ — 0, A± — .2,Bi = 
.3, mi = 2,m = 3,v t = 0), string 2 (a 2 = 2ir/3,A 2 = .l,B 2 = .2,m 2 = 3,ra 2 = 4,v t = 0) 
and string 3 (a 3 = 4n/3,A 3 = .2,B 3 = .4, m 3 = l,n 3 = 3,v t = 0). A snapshot of this loop 
shortly after the beginning of the simulation is shown in Figure El One can see the six kinks 
propagating away from the junctions. 
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Figure 7: A[a'- =2 ](z) as a function of z. We evolved until tf = 50 which corresponds to 
z f,2 = 47.4. Each point above zero corresponds to a left moving kink propagating on string 
number 2. As the system evolves, many kinks are created and although the vast majority has 
a very small amplitude, a large number has an amplitude of the order of 1. 

As one can see in Fig. [7] the total number of (left moving kinks) that have propagated on 
string number 2 is very large, and even though the amplitude of many of those kinks is small 
there still is a large number of those kinks with an amplitude larger than e.g. 1/4. This is in 



17 



line with the results in the previous section. We can be more precise by calculating Q ■ , the 
number of (left moving) kinks propagating on string j with an amplitude larger than 1/4, as a 
function of time measured in units L defined in ( 13"B"1) . The result is shown in Fig. El One clearly 
sees that after some initial fluctuations, a regime of exponential growth sets in. 

Remarkably, the slope of the curves in Fig. [S] is larger than the slope of the corresponding 
function logQ > shown in Fig. O for the simple model of a loop with junctions discussed in 
Section 4 (where we do do not take in account the loop evolution). The reason that loop dy- 
namics has this effect on the proliferation process is simply because the loop dynamics generally 
implies that one of the strings shrinks. On this string, kinks propagate more frequently back 
and forward between the junctions, thereby increasing the rate of proliferation. One expects 
therefore that, for a given set of tensions, proliferation is in fact the least efficient when all 
lengths are constant and equal, which is just the case considered in Section 4, and that the rate 
obtained here is more realistic. We note when the string that is shrinking becomes small, kinks 
tend to have a smaller amplitude on it. 
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Figure 8: The logarithm of the number of (left moving) kinks of amplitude A > 1/4 on each 
string, as a function of time measured in units of the total invariant length of the loop L defined 
in ( 138B (which remains constant in time), for the case of strings with equal tensions. After some 
initial fluctuations, a regime of exponential growth sets in. 



Finally, we note that these simulations end at tf = 50 because the junctions collide briefly 
after that time, as is evident from figure M At that time, the total number of kinks on this 
loop is of the order 10 4 . 
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Figure 9: Evolution of the "Invariant lengths" of the strings. The simulations end around 
tf = 50 when string 1 shrinks to a point and the junctions collide. 

5.3 Discussion 

The details of the evolution of the number of kinks evidently depend somewhat on the initial 
conditions. However, our central result that the number of large amplitude kinks proliferates 
exponentially appears to be universal. In particular, it is a robust feature of the evolution for 
the wide range of parameter values we have scanned. 

The loop evolution has two important implications for the evolution of kinks. First, as 
discussed above, it enhances the proliferation rate because the length of one of the strings 
generally decreases. On this string, kinks propagate more frequently back and forward between 
the junctions. 

However, at the same time this possibly provides an end to the proliferation process. Indeed, 
for all the initial conditions that we have tried, one of the strings always ended up shrinking to 
a point, resulting in the collision of the two junctions. The result of such a collision is unclear 
and depends on the physics of the underlying theory [25] • The junctions might disappear thus 
ending the proliferation, or two new junctions might form and proliferation may then continue. 

Gravitational backreaction (e.g. by rounding off the kinks or by inducing the decay of the 
loop due to important radiation from the many kinks) may also end the proliferation. Finally, 
radiation of other fields might also become important as the number of kinks increases and 
could also play a limiting role. 

6 Conclusion 

Motivated by their effects on the gravitational wave emission of cosmic strings and superstrings, 
we have studied the dynamics of kinks on strings with junctions. We have concentrated in 
particular on the evolution of the number of very sharp - or equivalently, large amplitude - 
kinks since it turns out these provide the dominant contribution to the GW burst signal from 
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kinks on a network of strings [IH] • 

The propagation of a kink through a Y-j unction leads to the formation of three daughter 
kinks - one reflected kink and two transmitted kinks. We first showed analytically that, for a 
specific initially static junction configuration, one or two of the daughter kinks can be sharper 
than the incoming kink. This turns out not to be an isolated case: the amplification of kinks 
through their interaction with junctions is a rather generic phenomenon. Indeed we showed 
in Section that, assuming a random distribution for the four incoming waves specifying the 
junction configuration and for equal string tensions, kinks are amplified in a significant region 
of configuration space. In appendix A we have generalized this calculation to strings of different 
tensions finding similar results. 

The dominant contribution to the GW signal from a network of strings comes from the 
loops. In Sections H] and \5\ we have therefore studied the evolution of kinks on loops with 
junctions. We have considered loops which, for simplicity, contain 2 junctions. If one neglects 
the loop dynamics and assumes that i) all strings joining the two junctions have essentially the 
same invariant length, and also that ii) when a kink of known amplitude reaches a junction, 
the configuration of the latter is randomly drawn among those that yield the correct amplitude 
for the original kink, then one finds the amplification rate is sufficient to sustain an exponential 
growth of the number of large amplitude kinks. For a wide range of tensions, the coefficient in 
the exponent appears to be of order one when time is measured in units of L. The origin of 
this exponential growth was illustrated with a toy model in Appendix B. 

We have included the effect of the loop dynamics on the proliferation and amplification of 
kinks in Section El where we numerically integrated the equations of motion of a loop with two 
junctions for a rather general class of initial conditions. Our simulations generalize those of 
[25] who considered an initially static and symmetric planar loop configuration. Interestingly, 
we find the loop evolution somewhat enhances the rate at which the number of large amplitude 
kink grows. To a large extent this difference can be traced to the fact that under evolution, 
generically one of the strings shrinks. On this string, kinks propagate more frequently back 
and forward between the junctions, thereby increasing the rate of proliferation and thus also 
the number of large amplitude kinks. By the time the junctions on the loop collide we typically 
find at least as many as ~ 10 4 large amplitude kinks in the equal tension case. 

We note, however, that our simulations do not take into account intercommutations between 
strings. These may lead to the creation of new junctions as well as several smaller loops (with 
many kinks). Further, we end our simulations just before the junctions collide and have not 
addressed the subsequent evolution of the system. Finally, we note that backreaction effects 
will become increasingly important as the number of kinks grows, and that our simulations do 
not include this effect. 

Nevertheless our findings suggest that if an evolving network of strings with junctions con- 
tains a population of loops with junctions, these typically contain a large number of very sharp 
kinks. The implications of this for the GW signal emitted by networks of this kind will be 
discussed elsewhere 
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Appendix A 

In this Appendix we study the dependence of the distributions p{Cj) of transmission coefficients 
on the ratios of the tensions of each of the strings, for three semi-infinite strings meeting at 
a junction. The case of equal tensions was discussed in the text in Section 3 where it was 
found that the distributions of both the reflected and the transmitted kinks have significant 
tails where Cj > 1. As in Sec 3, we consider a kink moving towards the junction on string 1 
and assume a uniform distribution on the unit sphere of the vectors b' that specify the junction 
configuration. 

Incoming kink on a light string: \i\ <C [12 ~ A^3 

The kink reaches the junction on a string which is much lighter than the other two strings. 
One might therefore expect, on average, the kinks transmitted to the heavier strings to have 
reduced amplitudes and the reflected kink to have an amplitude comparable or even enhanced 
relative to the amplitude of the incoming kink. 

This is indeed what we find. Figure HU] shows the marginal distributions p{Cj) for strings 
of tensions p,\ = 0.1, /i 2 = /i 3 = 1. One sees the distribution of C\ is relatively flat, with a 
probability of amplification P{C\ > 1) larger than ten per cent. On the other hand p(0 2 ) is 
sharply peaked around a value much smaller than one. 

The existence of this peak can be understood from eq. ( 1771) . Let e = fXi/fi <C 1 and take 
/i 2 ~ p,$ so that V\j ii ~ 0(e). Then from (jHJ) it follows that sf + s^ = 0(e) so that 
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Now generically one has 1 — s 2 ~ 0(e) for the above configuration of tensions [26j. Therefore, 
since the functions b' 2 and b' 3 do not change when the kink crosses the junctions, it follows that 

^[a 2 ] = ^||a 2 --a 2 + ||~0(e)A[b;] => 2 ~ 0(e). (43) 

The coefficient 2 will be large only if either si or s 2 is close to minus one. However the 
probability for this is small; according to [2S] only three per cent of the s 2 lie in the interval 
[—1,-0.9], which explains why very few events are seen in the tail at large values of the 
distributions of 2 and O3. 

Finally, the average values of the transmission coefficients for the above set of tensions are 
given by 

(Ci) = 0.49, (C 2 ) = (0 3 ) = 0.09. (44) 
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Figure 10: Distributions of the different 
transmission coefficients for tensions /ii = 0.1, 
^2 = A*3 = 1) with the incoming kink on string 
1. The average values are (C\) = 0.49, 
(C 2 ) = (C 3 ) = 0.09 and P(d < 1) = 0.88, 
P(C 2 < 1) = P(C 3 < 1) = 0.99. 
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Figure 11: The same distributions as in FigfTUI 
but for tensions [i\ = 1, /i 2 = 0.1 and // 3 = 1. 
Here (d) = 0.09, (C 2 ) = 0.72, (C 3 ) = 0.99 
and P(Ci < 1) = 0.99, P(C 2 < 1) = 0.81, 
P(C 3 < 1) = 0.53. 



Incoming kink on a heavy string: Hi ~ fj,^ ^> /i2 

In this case, the kink propagates on a heavy string towards a junction consisting of another 
heavy string as well as a very light one. The light string is expected to play a minor role while 
the heavy strings almost behave as a single long string without junction on which the kink 
simply propagates. Hence the reflected kink amplitude C\ should be peaked about a very small 
value, while the transmitted kink amplitude on the heavy string C% should be peaked about 1. 
This is indeed what we find, as shown in Fig fTTl) . and an analytic argument similar to the one 
given above can be used to explain the position of the peaks in the distributions^ of C\ and 
C 2 . 

We now have s\ + S3 = O(e). The argument for C\ is exactly the same as above, since 




2/i 3 1 - s 3 



O(e) 



Efc^fc l + si 

l+0(£) l+0(e) 



h> 



(45) 



For C3, we have 



S3 



2/13 



S3 



Efc/^fc, 



K- 



0(e) 




2/ii 



Si 



l+0(e) l+0(e) 



b; 



(46) 



The zeroth order in e does not vanish in Aa 3 , because b[ undergoes a jump. Instead we have, 
in generic configurations where the other coefficients are of order 1, Aa 3 rj Ab^ which means 
C 3 «l. 



lx The amplitude of the kink transmitted to the light string depends crucially on the configuration, as indicated 
by the relatively flat distribution of C^. 
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The mean values (Cj) for the set of tensions of Fig [TT] are smaller than 1; 
(d) = 0.09, (C 2 ) = 0.72, (C 3 ) = 0.99, 
but the probability of having an amplification on strings 2 and 3 is significant 
p(C t > 1) = 0.01, P(C 2 > 1) = 0.19, P(C 3 > 1) = 0.47. 



(47) 



(48) 



Again we note that, for string 3, although amplification is frequent this is mostly limited in 
amplitude since the distribution is sharply peaked around a value close to 1. 

Incoming kink on a heavy string with fix < ^2 + [13 

In the example shown in Fig[T21 fi 2 = Hz = 1 and /xi = 1.9. The (superimposed) distributions 
p(C 2 ), p{C$) are sharply peaked around 1, so the amplitude of the kinks transmitted to the 
light strings is comparable to that of the incoming kink on the heavy string. The distribution 
of C\ is much flatter, with a significant tail at large values. For this set of tensions one has 



P(d > 1) =0.11. 



(49) 



The presence and position of the peak in p(C 2 ) can again be explained using an analytic 
argument. Let e = - < 1. For most configurations, one has s 1 = — l+0(e), while s 2 = 1— O(e) 



and sf = 1 - O(e) 



(see also (jHD). Thus from eq. ( 12T1) . 



/± 




V2 \ fa/ 2^3 1 



'V.- ** ' 



/' 



s 2 



0(e) 



O(e) 




(50) 



Hence to zeroth order in e, Aa' 2 ?a Ab' x and therefore C 2 ~ 1. A similar argument applies to 

o 3 . 




Figure 12: Distributions of the different trans- 
mission coefficients for tensions fii = 1.9 and 
fi 2 = 1 = /X3. The mean values are (Ci) = 
0.49, (C 2 ) = 0.99, (0 3 ) = 0.99 and P{d < 
1) = 0.88, P{C 2 < 1) = 0.53, P(C 3 < 1) = 
0.53 




Figure 13: The same distributions for tensions 
/ii = 1 = /i3 and /i 2 = 1.9. Here (Ci) = 0.09, 
(0 2 ) = 0.72, (0 3 ) = 0.99 and P(Ci < 1) = 
0.99, P{C 2 < 1) = 0.81, P(C 3 < 1) = 0.99 
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Incoming kink on light string of tension \i\ > \x^ ~ A^3 

Finally in Fig[T3J we show the distributions of the transmission coefficients for the following set 
of tensions, \x\ = \x% = 1 and \x 2 = 1.9. One sees p{C\) and p{C%) are peaked at small values. 
The typical amplitude of the kinks reflected on the light string fxi and transmitted on the other 
light string of tebnsion /i 3 is very small. The distribution of C 2 is flatter, so the amplitude of the 
kink transmitted to the heavy string depends strongly on the configuration. This distribution 
also has a tail for C 2 > 1, so kinks are amplified in a substantial volume of configuration space: 

P(C 2 > 1) =0.19. (51) 

Again the presence and the position of the peak in p(Ci) and p(C 3 ) can be explained 
analytically. The argument is as above though now e = - <C 1 and in most configurations, 
si = 1 - 0(e) and s 3 = 1 - 0(e) while s 2 = -1 + 0(e). 

Appendix B 

Here we illustrate with a very simple toy model that an exponential behaviour of the kind 
exhibited in Sections H] and El is actually very generic. 

Consider the situation in which one constructs a tree, in which each node has two daughter 
nodes so that the nth generation therefore contains 2 n nodes. Each node contains an amplitude 
and the amplitude of each of its daughter nodes is obtained by multiplying this amplitude 
by a factor drawn randomly using a Bernoulli distribution: 2 with probability p or 1/2 with 
probability 1— p. This means that there is either an amplification (by a factor of 2) or a reduction 
of amplitude (by a factor of 1/2). Since we are interested in the case where amplifications are 
the least probable outcome, we set 

p < 1/2. (52) 

In order to initialize the experiment, we need to define the amplitude of the initial node (0-th 
generation): we set it to be 1. 




A A A 

Figure 14: Random experiment described in this section 

Let p n be the probability that any given node of the nth generation (say ki^ n for instance) 
has an amplitude larger than 1. Since p < 1/2, we expect p n to be small when n becomes large. 
Our goal here is to compute analytically (in the large n limit) p n . 

In order to be larger than 1, k\ yU must be the result of more amplifications than reductions 
so we clearly have 

Pn= E f?V(l-p) B -*- (53) 

k>n/2 V s / 
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Using Stirling's formula and transforming the sum into an integral using x = k/n, we can easily 
obtain 

p n = ^E f 1 X e nMx) dx (54) 

Jl / 2 yjx{\ ~ X) 

with 

f p (x) = -xln(-) - (1 - x) ln(^— £ )- (55) 

p 1 — p 

In the interval [1/2,1], f p (x) is a decreasing function (because p < 1/2) and therefore, it is 
maximal for x = 1/2. Because of the exponential in the integral, we then expect the dominant 
contribution to p n to come from the vicinity of 1/2. More precisely, by computing carefully the 
integral, it is possible to show that 



As expected, p n decreases exponentially since 2Jp(l — p) < 1. However, we are interested 
in the number N n of nodes of the n-th generation that have an amplitude larger than 1. If n is 
large enough, this number will typically be 2 n p n : 



N n ~ ,i „, 7= W-P) • (57) 



This number increases exponentially when AJp(l — p) > 1 i.e. p(l — p) > i. This condition is 

satisfied as soon as p > | — ^ ~ 0.07 (remember that we also imposed p < 1/2). 

This means that even if amplifications are not the most probable outcome of amplitude 
transmissions (as is the case for our physical system), the number of nodes of the nth generation 
that have an amplitude larger than some fixed value increases exponentially provided that p is 
not too small, i.e. provided that amplifications are not too rare. 

Note that the fraction of nodes that have an amplitude larger than 1 is given by p n and 
tends to exponentially fast and that the proliferation of large amplitude nodes is only possible 
because the total number (2 n in our example) increases faster than p n decreases. As expected, 
the very large majority of nodes have small amplitude. 

Of course, this very simple example is not a good physical picture of our kink proliferation 
(amplitudes can be larger than 1 in this example). Indeed, the Bernoulli distribution is a crude 
simplification of p^Ci, C2, C3) which does not take into account the fact that the amplitude 
of the daughter kinks are actually correlated and forgets about the amplitude of the incoming 
kink Ai. However, we expect the general mechanism to remain the same. 
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